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ABSTRACT 

We have analyzed a number of intense X-ray flares observed in the Chandra Orion Ultradeep Project 
(COUP), a 13 days observation of the Orion Nebula Cluster (ONC), concentrating on the events with 
the highest statistics (in terms of photon flux and event duration). Analysis of the flare decay allows 
to determine the physical parameters of the flaring structure, in particular its size and (using the peak 
temperature and emission measure of the event) the peak density, pressure and minimum confining 
magnetic field. A total of 32 events, representing the most powerful ~ 1% of COUP flares, have 
sufficient statistics and are sufficiently well resolved to grant a detailed analysis. A broad range of 
decay times are present in the sample of flares, with t\ c (the 1/e decay time) ranging from 10 to 
400 ks. Peak flare temperatures are often very high, with half of the flares in the sample showing 
temperatures in excess of 100 MK. Significant sustained heating is present in the majority of the 
flares. The magnetic structures which are found, from the analysis of the flare's decay, to confine 
the plasma are in a number of cases very long, with semi-lengths up to ~ 10 12 cm, implying the 
presence of magnetic fields of hundreds of G (necessary to confine the hot flaring plasma) extending 
to comparable distance from the stellar photosphere. These very large sizes for the flaring structures 
(L > J?,) are not found in more evolved stars, where, almost invariably, the same type of analysis 
results in structures with L < i?*. As the majority of young stars in the ONC are surrounded by disks, 
we speculate that the large magnetic structures which confine the flaring plasma are actually the same 
type of structures which channel the plasma in the magnetospheric accretion paradigm, connecting 
the star's photosphere with the accretion disk. 

Subject headings: magnetic fields - open clusters and associations: Individual (Orion) - planetary 
systems: protoplanetary disks -stars: activity - stars: pre-main sequence - X-rays: 
stars 



1. INTRODUCTION 

Beginning in their protostellar (Class I) phase, young 
stellar objects (YSOs) hav e long been known as copious 
sources of X-ray emiss ion ijFeigelson fc Montmerleiri999t 
iFavata fc Mice!all2003T) . The presence of frequent flares, 
observed both in accreting and non-accreting sources, 
shows that the emitting plasma is magnetically con- 
Electronic address: Fabio.Favata@rssd.esa.int 



fined, and in most cases also magnetically heated. In 
more evolved (Class III, non-accreting) YSOs the X-ray 
emission mechanism is likely to be very similar to the 
one operating in similarly active older stars, involving a 
scaled up version of the solar corona, with the emitting 
plasma being entirely confined in magnetic loop struc- 
tures and heated by similar mechanisms as in less ac- 
tive stars (most likely magnetic energy dissipa tion caused 
by the shear at the loop's footpoints, e.g. iPeter et all 
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12001 . In accreting sources (Class I and Class II), on the 
other hand, recent evidence points to the X-ray emis- 
sion mechanism b eing, at least in part, different: statisti- 
cally, in the ONC l|Flaccomio et al.l2003t iPreibisch et all 
2005) accreting YSOs are observed to be less X-ray lumi- 
nous than non-accreting ones, although the actual mech- 
anism causing this difference is not understood. Two 
accreting YSOs have been observed to date in X-rays at 
high spectral resolution: TW Hya, lIKastner et alJl2002t 
i Stelzer fe Schmittl l2004|) . and BP Tau (jSchmitt et all 
2005) . In both cases the O vn He- like triplet shows very 
low f/i ratios, which can be caused by very high den- 
sities and/or by a very high UV ambient flux. Such 
low f/i is not observed in any "normal" (ZAMS, MS or 
more e volved) coronal source. As e.g. discussed by Drake 
(2005) both effects are most likely associated with the 
accretion stream (obviously not present in non-accreting 
stars). While the X-ray emission from TW Hya is very 
soft (T ~ 3 MK), and thus could be explained as (largely) 
driven by the accretion shock, both cool and hot plasma 
(up to 30 MK) is present in BP Tau; the hot plasma can- 
not originate in the accretion shock (not enough grav- 
itational energy is available), and thus in these cases 
accretion-driven X-ray emission likely coexists with a 
more "canonical" corona. 

An obvious question to ask, also in the light of the re- 
cent results concerning TW Hya and BP Tau, is therefore 
whether the different types of YSOs have different types 
of X-ray emitting structures (and thus of confining mag- 
netic fields) in terms of size, location, structuring and 
density. One of the few diagnostic tools which can pro- 
vide an insight in this question is the analysis of intense 
X-ray flares. In this context, a flare is defined as a sudden 
impulsive rise in the plasma temperature, immediately 
followed by a rise in the X-ray luminosity of the source 
(typically by a factor of at least a few times the quies- 
cent luminosity), and later followed by a slower, roughly 
exponential decay of both the X-ray luminosity and tem- 
perature. These types of events a re often observed both 
in the Sun (see e.g. the review bv lPriest fc ForbeHl2002ft 
and in most active stars, and their study has been the 
subject of a copious litera ture (see e.g. the relevan t chap- 
ter reviewing the topic in lFavata fc MiceTall2003|) . 

The Chandra Orion Ultradeep Project (COUP) is an 
unique long (13 days span with 9.7 days of effective ex- 
posure) X-ray observation of the Orion Nebula Cluster 
(ONC). T he details of the CO UP observation are de- 
scribed bv iCetman et al.l (j2£)05). With its long, nearly 
uninterrupted observation, the COUP sample offers an 
unique opportunity to study flares in YSOs, and thus to 
study the type of magnetic structures present in young 
stars, and in particular in stars which are either actively 
accreting or which are surrounded by inactive disks. The 
magnetic field structure in YSOs, is expected to be differ- 
ent from the one in older, isolated stars. Magnetic field 
lines (and thus possibly magnetic loops) connecting the 
stellar photosphere and the disk have been postulated 
in th e context of the magnetospheric l|Hartmann et all 
1994) model of accretion, in which the accreting plasma 
originates from the disk and is channeled onto the stellar 
photosphere by the magnetic field. Magnetically fun- 
neled accretion both expla ins the observed broad op- 
tical emission line profiles (Ha rtmann et all 119941) and 
the ejection of high-velocity bipolar outflows l)Shu et aTl 



2000). 

However, no direct observational evidence for such 
extended magnetic structures has been available up 
to now. Strong magnetic fields have been measured 
in the photosphere of YSOs using Zecman splitting, 
and fields strengths up to 6 kG have been deter mined 
(see the review of IJohns-Krull fc Valentil |005; also 
IJohns-Krulletal] l2004TT ~ Such measurements cannot 
however determine the field's structure, and whether it 
is dominated by a large-scale, dipole structure, or by 
higher-order multipoles imposing a smaller scale to the 
field. In particular, the measurement of photospheric 
fields does not allow determination of how the field ex- 
tends into the circumstellar region, and in particular into 
the disk. Magnetic flux tubes are required by the mag- 
netospheric model of accretion to channel accreting ma- 
terial from the inner edge of the disk to the photosphere. 
Analysis of the flares' decay can help to determine if the 
same flux tubes are also the seat of more energetic pro- 
cesses, i.e. if the plasma is heated to temperatures re- 
sulting in X-ray emission and in the attendant flaring 
activity. 

A number of intense X-ray flares on PMS stars have 
been previously observed and studied. Prior to the 
XMM-Newton and Chandra observations some 6 flares 
on PMS stars had been a nalyzed in detail (see review 
by iFavata fc Micelal I2003D , although a large number of 
fainter events (for which therefore physical parameters 
could not be derive d) had also been observed - see e.g. 
IStelzer et all (2000) for a survey of flares in th e Taurus- 
Auriga-Perseus complex. IFavata et al.l l)2001f) analyzed 
four of the bright events on YSOs known at the time, 
finding, in all cases, that the flaring loop structures were 
of modest size (L < i?*), similar to what observed in 
active main sequence and more evolved stars. 

A number of flaring events on YSOs have been 
observed with Cha ndra and XMM-Newton (see e.g. 
Ilmanishi et al.ll2003l for a Chandra study of flares in the 
p Oph region ) , but very few have been analyzed in detail. 
iCrosso et al.l l)2004f) analyzed a flare observed with Chan- 
dra on the YSO LkHa312 in M78, also finding modest- 
sized loop structures (L — 0.2-0.5 R*). 

Therefore, the analysis of the flares observed up to now 
on YSOs has resulted in magnetic structures which are 
in all respects similar to the ones found in older stars, 
supporting a view in which the coronae of YSOs are sim- 
ilar to the ones of their more evolved counterparts, and 
supplying no evidence for the long magnetic structures 
expected to connect the star with its disk. As it will be 
shown in this paper, the analysis of the large flares occur- 
ring in the COUP sample provides a different picture: in 
addition to flares occurring in star-sized or smaller struc- 
tures, similar to what found up to now in both PMS and 
older stars, a number of events are found taking place 
in very long magnetic structures, with lengths of several 
stellar radii. As discussed in Sect. HO loops of this size 
around fast rotating stars are unlikely to be stable if an- 
chored on the star alone, and therefore such structures 
are likely to be extending from the star to the disk, pro- 
viding observational evidence for the type of structures 
postulated by magnetospheric accretion models. 

While COUP is the best sample available to date for 
the study of the size of coronal structures on YSOs, its 
limitations should also be stressed. The main one is the 
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relatively limited photon statistics of Chandra observa- 
tions of YSOs in the ONC: only the most intense and 
longer lasting flares will have sufficient statistics to al- 
low for the time-resolved spectral analysis needed by the 
approach used here. Therefore, the physical parameters 
derived in the present paper are not necessarily represen- 
tative of the "average" flaring structure in the coronae 
of YSOs (in the ONC or elsewhere); rather they must 
be interpreted as the structures associated with the most 
intense flares present in these coronae. Structures associ- 
ated with less intense flares (most likely smaller magnetic 
loops than the ones associated with very intense flares) 
are certainly present, and given the large large number 
of shorter and less intense flares visible in the light curve 
of these stars (which cannot be analyzed in detail due 
to the low statistics) these are common. Indeed, even 
in our biased sample a number of compact flaring events 
is present, and, like what observed in older stars, most 
likely confined in loops anchored on the stellar photo- 
sphere, with no interaction with the disk. 

The present paper is structured as follows: the selec- 
tion of our sample is described in Sect. |5J the analysis 
procedures adopted in Sect. EH the results from the analy- 
sis are presented in Sect.0J which includes detailed mod- 
eling of a sample flare (on COUP 1343, Sect. I4.1|l . An 
example of a complex flaring event (COUP 450, which 
cannot be analyzed with the method adopted here) is dis- 
cussed in Sect. and the implications from the present 
work are discussed in Sect. El The conclusions are sum- 
marized in Sect.0 while Appendix [3 discusses in detail 
a number of notable individual events. 

2. SAMPLE SELECTION 

We have studied all flare s in the COUP samp le (de- 
scribed in detail elsewhere, iGetman et al J 1200 5 ) which 
have sufficient statistics for the analysis performed here. 
From a set of simulations we initially determined the 
minimum number of photons in a spectrum necessary to 
derive reliable spectral parameters given the typical X- 
ray spectra of COUP stars. A minimum of 750 photons 
was found to represent a reasonable compromise between 
the need for statistics and being able to analyze as large 
as possible a number of events. 

To subdivide the light curve in different interval, a 
Maximum Likelihood (ML) a lgorithm was used, the same 
as also used and discussed by jWolk et alj |2Pj05j[). Briefly, 
the rationale behind the ML approach is to recursively 
subdivide the light curve in a number of segments, with 
the property that the source's count rate is compatible 
with being constant during each segment. The minimum 
number of photons which are comprised in a given ML 
block is a configurable parameter. While for many appli- 
cations allowing blocks with as little as 1 photon is useful 
(e.g. to allow to track fast variability), for our analysis 
a minimum of 750 photons per ML blocks has been im- 
posed (to allow for meaningful spectral analysis). In all 
the following whenever the term "ML block" is used, it 
implies ML blocks with a minimum of 750 photons each. 

The ML block approach identifies (as discussed by 
iWolk et alJ l2005|) a 'characteristic' level, essentially a 
quiescent level during the observation, as well as blocks 
which are significantly above this level and during which 
the X-ray luminosity is constant. Note that our con- 
straint of > 750 photons per block results, in a number 
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Fig. 1. — Top panel: the COUP light curve of source 1343; 
the ML blocks at the characteristic level are plotted in light blue 
and the ML blocks with a higher level are plotted in orange. The 
source's ACIS count rate is given in the left axis, while the right axis 
gives the approximate number of counts in each 2 ks bin. Time 
is given in seconds from the beginning of the observation in the 
bottom axis, and in calendar days (January 2003) in the top axis. 
Bottom panel: evolution of the flare. The ML blocked light curve 
with only the flaring event shown is plotted in the bottom right 
quadrant, while the temperature and emission measure evolution 
are plotted in the bottom left and top right quadrants respectively. 
The top left quadrant shows the evolution of the flare in the the 
log T— log V EM plane, with the points connected by a das hed line 
to allow to follow the event's evolution. See the text (Sect. HOI for 
a complete description of the diagrams. 

of cases, in significant "under-binning" of the light curve, 
i.e. the light curve is clearly variable within a given block. 
This is for example visible in the case of COUP 649 (see 
Sect. EH) - 

Candidate flares for analysis have been initially hand- 
selected by inspection of the atlas of light curves pro- 
duced by the ML blocking for the complete COUP sam- 
ple, and choosing all sources which had at least three con- 
tiguous ML blocks above the characteristic values. These 
are candidate flaring events with sufficient information 
to grant a more detailed analysis. For example the light 
curve of COUP 1343 is plotted in Fig.^ where (as in all 
other light curves) the ML blocks at the characteristic 
level are plotted in light blue and the ML blocks with a 
higher level are plotted in orange. The latter ML blocks 
form the candidate flaring events. For COUP 1343, in 
addition to the flare being analyzed (at the beginning of 
the observation) there are additional ML blocks above 
the characteristic level, which have been ignored as they 
are not part of the main flare. This is also common in a 
number of other sources. 

Following the spectral analysis described below, the 
evolution of each candidate flaring event in the logT- 
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log V EM plane (a proxy for the logT-logy 7 ?^ plane, as 
described in Sect.OJ has been examin ed. Events evolving 
according to the pattern discussed bv lReale et al.l l)1997fl 
show, after the flare peak, a decay of both the tempera- 
ture and the emission measure, more or less regular de- 
pending on the event characteristics. The event shown 
in Fig^shows this 'canonical' behavior: as visible in the 
the logT vs. log \J EM diagram (where the event evolves 
temporally in the clockwise direction), the temperature 
peaks while the emission measure is still rising (block 
2), and has already started to decay when the emission 
measure peaks (block 3). Then (blocks 4-9) both the 
temperature and the emission measure decay regularly. 
Two smaller flares are also visible in the light curve, one 
on day 13 and the other on day 21. As not enough statis- 
tics is available in each case, they have not been further 
considered. 

We discarded for the purpose of the present analysis 
all events for which the evolution did not proceed in a 
sufficiently regular way (e.g. events for which the tem- 
perature did not decay). These can either be flares with 
different underlying mechanisms, which cannot therefore 
be analyzed within the framework used here (such as e.g. 
the two-ribbon flares seen in the Sun) or different types of 
variability (e.g. induced by rotational modulation, or by 
the emergence of new active regions, etc.). Only events 
which could be analyzed in an useful way (i.e. produc- 
ing well constrained physical parameters within the as- 
sumed framework) have been retained in the final sample. 
One example of a particularly complex and intense flare, 
which cannot be analyzed with the approach used here, 
is shown and discussed in Sect. 03 

Clearly, as mentioned in Sect. ^ such procedure does 
not produce an unbiased sample. The most obvious bias 
is toward events which have sufficient X-ray flux and du- 
ration to yield enough time-resolved spectra to allow the 
analysis. The resulting physical parameters are there- 
fore unlikely to be representative of the 'typical' flaring 
structures present in YSOs in the ONC, but will rather 
represent the more extreme cases (and thus the largest 
magnetic structures present). Quantifying the biases in- 
troduced by our selection procedure is not obvious. An 
indicati on can however co me from comparison with the 
study of lWolk et all 1120051) . who have performed an anal- 
ysis of flares in young solar analogs in COUP using an 
automatic procedure. The young solar analog sample 
consists of 28 out of 1616 COUP sources. A total of 
41 distinct flares were identified on 26 of them, for an 
average of ~ 1.5 flares per stars during the COUP ob- 
servation. If the statistics on solar analogs were valid for 
the complete COUP sample (an assumption which needs 
to be verified) this would imply some 2400 flares in the 
observation. The 32 events analyzed here represent only 
1% of the COUP population of flares expected extrapo- 
lating the statistics for solar-mass stars. We are therefore 
studying (by necessity, imposed by photon statistics) the 
tip of an much larger iceberg. 

3. ANALYSIS METHOD 

The analysis of the decay of flares is a classic tool to 
derive the size of the flaring structure, and thus by infer- 
ence other quantities such as the plasma density and the 
confining magnetic f ield (as e.g. discussed for the Sun by 
iCulhane et al.lll970]) . Stellar flares have been observed 



since the beginning of stellar X-ray astronomy, and they 
have been analyzed in a variety of ways, most of which 
proceeded by analogy to the characteristics of (spatially 
resolved) solar flares. In the case of the Sun, flares occur 
in localized regions in the corona, often in single magnetic 
loop structures, which can remain mostly unchanged dur- 
ing the flare evolution. 

While the heating mechanism causing the flare still 
remains not understood in its details, it is generally as- 
sumed to be unlocked by the release of energy stored in 
the magnetic field. The impulsive energy release gen- 
erates either a strong thermal front or particle beams, 
which, channeled by the magnetic field, hit the photo- 
sphere and cause plasma evaporation; the plasma fills 
the coronal loop(s), and produces the observed increase 
in X-ray luminosity. The decay of the temperature and 
the X-ray luminosity is caused jointly by the radiative 
losses and by heat conduction back to the photosphere. 

Both cooling mechanisms depend on the geometry of 
the flaring structure, and thus the decay time scale of 
flares is related to the size of the coronal structure (or 
structures). Stellar flares have been used to infer coronal 
structure size s since the early stages of coronal X-ray as- 
tronomy (e.g. iPallavicini et al-l fl 990) . However, up and 
until the late '90s, it was largely assumed that the heating 
event was impulsive, and concentrated at the beginning 
of the flare, coincident with the observed impulsive rise in 
temperature, so that long decay times always led to infer 
large coronal structures, in the form of very long loops. 
In the '90s studies of solar flares showed that the heating 
mechanisms often extends well into the decay phase, so 
that sustained heating and thermodynamic cooling pro- 
cesses are competing, in some cases resulting in long de- 
cay times also for flares confined in compact structures. 
Thus, analyses of stellar flares based on the assumption 
of impulsive heating would likely result in longer coronal 
structures than actually pres ent on the star (see fo r a 
d etailed review of past r esults iFavata fc Mice' la 2003). 

iSvlwester et all (I1993T) showed that, in solar flares, the 
slope of the flare decay in the temperature-density plane 
is a se nsitive diagnostic o f the presence of sustained heat- 
ing. iReale etall <|1997|) later extended the approach 
showing that the slope ( of the flare decay in a logT vs. 
log n e diagram gives a quantitative measure of the time- 
scale of sustained heating, which can be used to correct 
the observed decay time scale and to derive the intrin- 
sic (thermodynamic) decay time, thus allowin g us to ob- 
tain t he actual size of the flaring structure. IReale et alJ 
1)1997(1 developed a methodology for the analysis of stellar 
flares, based on grids of hydrodynamic models of flaring 
loops, which was verified on solar flares (for which the 
size of individual flaring structures can be obtained from 
imaging observations), and which has been applied to a 
significant number of stellar events on a variety of stars, 
including YSOs. In the vast majority of cases, the anal- 
ysis of stellar events has shown that sustained heating is 
present, so that flaring structures in stellar coronae are 
generally smaller than previously thought. 

3.1. Uniform cooling loop modeling 

The rise time in both solar and stellar flares is almost 
invariably observed to be significantly shorter than the 
decay time, and therefore it has been often assumed, in 
deriving the flare's physical parameters, that the heat- 
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ing is impulsive, concentrated at the beginning of the 
event, and that the decay takes place in an "undisturbed" 
fashion from a loop in near-equilibrium at the flare's 
peak, according to the thermodynamic cooling time of 
the plasma. The two cooling processes which determine 
the decay time of the flare are thermal conduction down- 
wards to the chromosphere and radiation, each with its 
characteristic 1/e decay time, 



Tcond 



3nkT 
kTV 2 /L 2 



7"rad 



3nkT 
n 2 P(T) 



(1) 



where n is the plasma density, k is the thermal conduc- 
tivity and P(T) is the plasma emissivity per unit emis- 
sion measure. Both decay times depend on the loop's 
length, T con d explicitly through the L 2 term, r ra( j im- 
plicitly through the density's dependence. The effective 
cooling time of the loop is the combination of the two, 
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More generally it has been shown bv lSerio et alJ j|1991j) 
that the decay time of a flaring loop starting from equi- 
librium and decaying freely is linearly related to its half 
length, through 

r th-\/Tpk 



3.7 x io- 



ta) 



where L is in units of 10 9 cm, T t h in sec and T p ^ is the 
peak temperature of the plasma in the flaring loop in 
units of 10 7 K. A fast decay thus implies a short loop, 
while a slow decay implies a long loop. In all the follow- 
ing, unless otherwise specified, the term "length" will be 
used to indicate the loop's half length, from one of the 
footpoints to the loop apex. 

In the 80's and 90's a number of approaches to the 
analysis of stellar flare decays have b een used (as dis- 
cusse d in detail by Reald l2002t see also lFavata fc Micelal 
l20nH for a review of previous literature). While differing 
in their details, they were mostly equivalent to the use 
of Eq. 13 above, and their key assumption was that the 
flaring loop always decays, after a short heating episode, 
in an undisturbed fashion. 

If heating does not switch off abruptly but rather de- 
cays slowly, the plasma is subject to prolonged heating 
which extends into the decay phase of the flare. As a 
result, the actual flare decay time determined from the 
flare's light curve, t\ c , will be longer than the intrinsic 
decay time r trl , and simple application of Eq. [21 will re- 
sult in an over-estimate of the size of the flaring loop. 
If sustained heating is present the observed flare decay 
time must be "corrected" to obtain a reliable estimate of 
the flaring region's size. 

In the case of solar flares lSvlwester et al.l l)1993|) showed 
that the slope £ of the flare decay in a logT vs. logro e di- 
agram provides a diagnostic of the presence of sustained 
heating, with a shallower slope (slower temperature de- 
cay) indicating an event with strongly sustained heating. 
This ap proach was e xtend ed to the analysis of stellar 
flares bv lReale et alJ (|1997f) , who showed that £ provides 
a quantitative diagnostic of the ratio between the intrin- 
sic and observed decay times, i.e. 



— = F(0 

Tth 



(4) 



The approach of iReale et al.l j|1997f) was based on nu- 
merical simulations of flaring loops, under the assump- 
tion of an exponentially decaying heating function (with 
a time scale th), and it was validated by comparing its 
results with the measured size of a number of (spatially 
resolved) solar flares. 

The actual functional form and parameters of F(() de- 
pends on the bandpass and spectral response of the in- 
strument used to observe the X-ray emission from the 
flare, and therefore needs to be separately determined 
for each instrument. Note that for stellar observations no 
density determination is normally available; in this case, 
the quantity V EM is used as a proxy to the density, un- 
der the assumption that the geometry of the flaring loop 
does not vary during the decay. 

For ACIS observations the relationship between £ and 
the ratio between the observed and intrinsic flare decay 
time has been calibrated for the present analysis, result- 
ing in 

? - no = + i.4i (5) 



nil 



£-0.32 



which is valid for 0.32 < £ < 1.5. The limits of applicabil- 
ity correspond on one side (£ ~ 1.5) to a freely decaying 
loop, with no heating (th = 0), to the other (£ = 0.32) to 
a sequence of quasi-static states for the loop (th 3> Tth), 
in which the heating time scale is so long as to mask the 
loop's intrinsic decay. Note that while it is not possible 
to have loops with £ > 1.5, flares with £ < 0.32 are pos- 
sible if the heating is not simply exponentially decaying 
but if it increases again after an initial decay. 

The temperature in Eq. [3] is the peak temperature 
in the flaring loop, which will contain plasma with 
a distribution of temperatures. The relationship be- 
tween the maximum temperature present in the loop 
and the temperature determined by performing a sim- 
ple 1-temperature fit to the integrated X-ray spectrum 
emitted by the loop has also been calibrated here. In the 
case of observations performed with ACIS is 



T pk = 0.068 x T r 



obs 



(0) 



where both temperatures are in K. 

Putting together the above elements, the length of the 
flaring structure is therefore given by 



L = 



3.7 x 10- 4 F(Q 



(7) 



Once the length of the flaring loop has been derived, 
it is possible to infer additional physical parameters of 
the flaring regions, with some additional assumptions. 
Given (from the analysis of the peak spectrum) a peak 
emission measure for the flare, a density can be derived 
if a volume is available. The above analysis only gives 
the loop length; if an assumption is made about the ratio 
between the loop's radius and its length (/3 = r/L) one 
can derive a volume, 



V = 2irf3 2 L 6 
and the resulting plasma density at flare peak 



EM 
~V~ 



< EM 
2irf3 2 L 3 



(8) 



(9) 
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The minimum magnetic field necessary to confine the 
flaring plasma can be simply estimated as 

B ~ v / 87rp = VSTrkriff (10) 

where p is the plasma pressure derived from the density 
and the temperature. 

For the solar corona (the only one accessible to imaging 
observations) typically (3 ~ 0.1, and this value has been 
used in most previous analyses of stellar flares. We have 
also assumed, in all of our estimates, (3 — 0.1, unless 
otherwise specified. The consequences of differences in 
j3, depe nding on the loops characteristics, are discussed 
in 

Sects. O and H 

3.2. Uncertainties and error analysis 

Both statistical and systematic uncertainties are of 
course present in the length estimate for the flaring loop 
obtained through Eq.[7] The main statistical uncertainty 
comes from the uncertainty in the determination of £ and 
therefore of F(£). Uncertainties in the determination of 
the peak temperature and of the decay time are usually 
significantly smaller than the uncertainty deriving from 
the .F(C), which, given the hyperbolic form of the func- 
tion, becomes potentially very significant at the lower 
end of the range of validity of £. In our analysis, the 
error bar reported for the length estimate is the one ob- 
tained by propagation of the statistical errors on £, T p k 
and ric, and it is dominated by the error on £. 

One possible additional source of uncertainty would 
be the lack of accurate knowledge about the star's grav- 
ity. Gravity may in principle bring systematic effects 
to the diagnostics, because it makes the plasma drain 
faster during the decay phase, flattening the flare's de- 
cay in the log T-log V EM plane, mimickin g the presence 
of sustained heating l|Reale fc Micelall998fl . In a number 
of cases, no estimates of either the radius or the mass of 
the parent star are known, so that no determination of 
the surface gravity is possible. However, given the low 
values of surface gravity in YSOs (with their large radii, 
typically g ~ 0.1 g 1 ©), the reduction of the gravity at dis- 
tances comparable to the stellar radius from the surface, 
and the very high temperatures observed at flare peaks, 
the pressure scale height of the flaring plasma is likely to 
be always larger than the length of the flaring structures, 
so that gravity is unlikely to be an additional source of 
uncertainty. 

More subtle is the issue of systematics, i.e. of uncer- 
tainties deriving from assumptions in the models which 
are not reflected by reality. The modeling of the flaring 
plasma structure is based on relatively simple physics, 
robust in its assumptions. The modeling is simplified 
by the one-dimensional nature of the problem, with the 
magnetic field confining the plasma except along the field 
lines, so that the precise shape of the flaring structure is 
not important. Another assumption is that the sustained 
heating is an exponentially decaying function; of course 
this is a simple parameterization which has no underly- 
ing physical reasoning behind it. However, as long as 
the temperature indeed decreases, and the slope in the 
log T-log \JEM plane is steeper than the one of the quasi- 
static sequence locus (C > 0.32) the heating must be a 
decreasing function of time. Its detailed functional form 
in this case has a limited impact on the loop length, so 



that the resulting values are robust agains t the details 
of the heating function. Also, as shown bv lReale et alJ 
(2002) strong, long-lasting flares must be strongly con- 
fined by the magnetic and would decay much faster than 
observed if the plasma were to break confinement. 

The derived parameters (density, magnetic fields) are 
subject to an additional uncertainty relative to the loop's 
aspect ratio (3. Our simple assumption about (3 clearly 
cannot hold for a loop with constant cylindrical cross 
section when this has a length of several stellar radii (as 
is the case in some of the events analyzed here) . In this 
case the footpoints would cover the whole stars or be 
larger than it, a clearly unrealistic configuration. Thus, 
such loops must, at the footpoints, have smaller aspect 
ratios. At the same time, though, a loop of such length is 
likely to follow a dipole-like magnetic field configuration 
at large distances from the star (with the higher-degree 
multipole terms becoming less important with distance) . 
This implies that it will likely have an expanded cross- 
section at the apex (where most of the flaring plasma will 
be) than at the footpoints. These two effects compensate 
each other (although to an unknown degree), and we have 
thus decided to still maintain (3 — 0.1 in the analysis, 
keeping in mind that for the longer loops the "effective 
/3" may be somewhat smaller. 

3.3. Application of the analysis to the COUP data 
To determine the evolution of COUP flares in the log T 
vs. log V EM diagram, the temperatures and emission 
measure of the flaring component have been separately 
determined for each ML block. All spectral fits were 
performed within the XSPEC package, using the mekal 
spectral emissivity model for coronal equilibrium plasma, 
for consis tency with the rest of the COUP analysis as de- 
scribed in lGetman et al.l l)2005f) . An example of the anal- 
ysis is shown in Fig.^ the top panel shows the raw ACIS 
light curve (at a resolution of 2 ks), with superimposed 
the light curve rebinned in ML blocks. The segments in 
light blue are the ones representative of the "character- 
istic" level for the source's emission, while the ones in 
orange are the ones significantly above the characteris- 
tic level, and which form the basis for the flare analysis. 
The count rate for the light curve is given in cts/ks on 
the left vertical axis, and in (approximate) counts per 2 
ks bin on the right vertical axis. The time scale is given 
in ks from the observation start on the bottom axis and 
in calendar days (January 2003) in the top axis. Also 
indicated are the radius and mass of the star (if known, 
in solar units) and the total number of photons in the 
light curve. 

In a number of cases the flaring source showed signif- 
icant pile-up in the ACIS detector during the flare. In 
this case the pile-up free light-curve was examined (ob- 
tained by extracting the photons from a r ing around the 
peak of the PSF, as described in Sect. 6 of lCetman et alJ 
2005). However, in most cases the statistics of the pile-up 
free light curve are so much lower that insufficient pho- 
tons were left for a detailed analysis of the flare. Only 
in the case of COUP 891 (Fig. Eft could the flare be an- 
alyzed, in spite of the presence of pile-up. The flares in 
COUP 107, 245,290, 394, 430 and 881, which would have 
had sufficient statistics and sampling, were lost from our 
sample because of pile-up. 

The source's "characteristic spectrum" has been deter- 
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mined by performing an absorbed multi-temperature fit 
to the integrated spectrum of all ML blocks compatible 
with the characteristic level of the source. The purpose of 
this fit was only to provide a noise-free, phenomenological 
description of the source's characteristic emission, which 
is effectively a form of background to be subtracted from 
the flaring emission. Our fit provides a model of the back- 
ground to be subtracted, which is more robust and less 
sensitive to statistical fluctuations than the subtraction 
of the characteristics spectrum itself, which in a num- 
ber of cases has rather low statistics. Therefore enough 
temperature components have been included until the 
fit has yielded a satisfactory x 2 , an d individual metal 
abundances have been left free to vary, independent of 
whether their best-fit value was physically meaningful or 
not. 

For each ML block with emission above the charac- 
teristic level, a spectral fit to the emission in excess of 
the characteristic level has been performed. In addi- 
tion to the characteristic component model with frozen 
fit parameters, we include a single additional absorbed 
thermal component, with global metal abundance set to 
Z = 0.3 Zq, thus determining the emission measure and 
temperature of the flaring component. The absorbing 
column density A(H) was left free to vary. Previous 
works on CCD resolution X-ray spectroscopy of very ac- 
tive stars, including PMS stars, shows that they typi- 
cally have low coronal abundance, around Z — 0.3 Zq, 
which we have therefore chosen here. While coronal 
abundance variations have been observed during intense 
stellar flares, at the temperatures of interest here the 
plasma is fully dominated by continuum emission, and 
therefore the exact value of the metallicity - which of- 
ten cannot be well constrained - has little influence on 
the best-fit temperature and emission measures. In all 
cases this simple onc-tcmpcraturc fit resulted in a satis- 
factory description of the spectrum, although a number 
of events have spectral peculiarities (regarding variations 
of abundance and absorbing column density) which will 
be discussed in detail in a future paper. These do not 
affect any of our conclusions. 

The bottom four-quadrant panel in Fig. ^ shows the 
results of our analysis in the case of COUP 1343. The 
bottom right quadrant shows the ML blocked light curve 
with the time axis expanded for clarity (only the flar- 
ing event is shown, and none of the characteristic level 
blocks), while the upper right and lower left quadrants 
show, respectively, the time evolution of the emission 
measure and of the temperature of the flaring compo- 
nent, as determined from the analysis of each of the ML 
flaring blocks, in a semi-log scale. The best fit exponen- 
tial to the decay is also plotted for both quantities as a 
continuous line. The top left quadrant shows the evolu- 
tion of the flare in the the log T-log \J EM plane, with 
the points connected by a dashed line to allow to follow 
the event's evolution. Points are numbered accordingly 
to the sequential number of the ML block. The blocks 
which have been used to determine the flare's evolution 
are marked by a dot, and the best-fitting decay is plotted 
as continuous line. Its slope C is the one used in Eq. [SJ 
For all such plots we have chosen to adopt the same scale 
in both axes, to allow an immediate comparison of dif- 
ferent events. 



In Fig. we plot four representative ACIS spectra for 
COUP 1343, illustrating the type and quality of the spec- 
tra used in the present paper. The top left panel of Fig. [21 
shows the characteristic spectrum of COUP 1343, with 
the 3-tempcrature fit used to describe it. This spectrum 
has been integrated for ~ 540 ks, and has 3821 source file 
counts. The other panels show (clockwise) the spectrum 
at the beginning of the flare's rise (block 2, when the 
plasma is hottest), at the beginning of the flare's decay 
(block 4) and at the end of the decay (block 9, when the 
plasma has already cooled to a temperature very similar 
to the quiescent value). The flare spectra contain 738 
photons over 4.7 ks, 1058 photons over 5.6 ks and 1214 
photons over 36 ks, respectively. 

For a large number of sources (14 out of 32), the peak 
temperature derived by the above procedure converges 
to a value higher than 10 keV, although invariably with 
large error bars. Detailed analysis of the fits shows that 
the fit procedure cannot determine with good accuracy 
the temperature of thermal components hotter than ~ 
100 MK on ACIS spectra; when hotter plasma is present, 
statistically equivalent fits can be obtained by arbitrarily 
fixing the temperature for this hot component to any 
value between 100 and 600 MK. We have thus decided to 
adopt a conservative approach, and for all cases for which 
the maximum best- fit temperature T b s > 100 MK the 
value 100 MK has been used in Eq. El (yielding T pk = 270 
MK). This will result in smaller loop sizes derived from 
Eq. than it would be the case if the actual maximum 
peak temperature were used. However the dependence 
of L on Tpk is relatively weak (going as ^/T p k); thus, 
even if the peak temperature were to be underestimated, 
in a few cases perhaps by up to a factor of two, this 
would result in an additional uncertainty in the size of 
the flaring loop of ~ 40%. 

4. RESULTS 

The physical parameters derived for each flare are 
listed in Table ^ F° r each source we report the max- 
imum observed temperature T \, s , the peak X-ray lumi- 
nosity Lx, the peak temperature at the loop apex T p k 
determined from Eq. H3 the slope £ of the flare decay 
in the logT vs. log V EM diagram, the ratio of the ob- 
served to the intrinsic decay times F(Q (determined from 
Eq. EJ, the observed decay time scale (in ks) t\ c deter- 
mined by fitting the flare decay binned in ML blocks. 
The length of the flaring loop L (from Eq. [7J is given 
in units of 10 10 cm, as well as, for stars which have a 
radius estimate published, in units of the stellar radius. 
The density n e computed from the peak emission mea- 
sure and from the volume of the flaring loop (Eq. EJ) is 
given in units of 10 10 cm -3 , and the minimum magnetic 
field necessary to confine the flaring plasma at the loop 
apex B (Eq. IT0|) is given in Gauss. Also given (when 
known) are the equivalent width of the Can IR triplet, a 
diagnostic of active accretion processes, and whether the 
source has A(K — L) and A(7 — A") excesses. Note that 
the peak X-ray luminosity is computed after subtraction 
of the characteristic emission from the flaring one. 

The derived quantities (n e ,B) are computed assum- 
ing (3 = 0.1. For long loops (of several stellar radii) 
this assumption may be incorrect, as such large loops 
would have footpoints covering much of the stellar sur- 
face (even though loop expansion may partially offset 
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TABLE 1 

Flare parameters from the ONC COUP sample 
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Note. — Sec main text (Scct.^ for a description of individual columns. The value for £ is given together with the uncertainty range 
deriving from the fit to the flare decay in the log T— log V EM plane. When £ is derived from only 2 decay points no uncertainty is given 
in the Table. In this case no uncertainty range is given for L; this value, as well as the derived values for n e and B : are in this case to be 
treated with caution. 

a C was derived from only two points in the flare decay and therefore no uncertainty range is given. 



this, see Sect. 13.21) . However, the dependence of n e and 
B on (3 is moderate, with n e oc /3" 1 and B oc f3~ ' 5 . 
Thus, the values in Tabled can be easily scaled to any 
value of /3. 

The flares analyzed here are all X-ray bright, but not 
exceptionally so, when compared to other intense flares 
already observed in YSOs. The peak X-ray luminosities 
reported in Tabled determined by fitting the spectrum 
of the peak ML block for the flare, range from 7 x 10 30 erg 
s _1 up to 8x 10 32 erg s _1 for the most X-ray bright event. 
The latter one takes place on COUP 1568 (discussed in 
detail in Appendix UJ), at M — 2.6 Mq one of the most 
massive stars in our sample. Thanks to the large off- 
set angle, the source does not suffer from pile-up even 
if the peak ACIS count rate approaches 1 cts/s. While 
certainly very high, the peak luminosity of 8 x 10 32 erg 
s _1 is comparable to the values observed in other PMS 
flares, e.g. th e 10 33 erg s -1 determ ined at peak in a sim- 
ilar band bv iTsuboi et al.l l|1998f) for the flare observed 
with the ASCA satellite on the weak line T Tau V773 
Tau. 

More remarkable are the high temperatures present in 
COUP flares. Up to recently, flaring events with tem- 
peratures reaching up to 100 MK where considered ex- 
ceptional. In the COUP sample, on the other hand, very 



high temperatures are common. While the ACIS spectral 
response does not allow to properly constrain tempera- 
tures in excess of 100 MK, the best fit values for the 
hottest ML block in the flare exceed 100 MK in about 
half of the events studied here. The fact that the tem- 
perature decay shows a clear regular pattern even when 
the best-fit values are above 100 MK (as for example is 
the case for COUP 1568), gives confidence in the fact 
that the peak observed temperatures are, in a number of 
cases, very high, hundreds of MK, even though their val- 
ues cannot be accurately determined using ACIS spectra 
(see Sect. ESJ). 

For the ONC YSO sources studied here, the nominal 
sizes of the flaring structures vary from a fraction of the 
stellar radius (e.g. COUP 7, 90, 141), with absolute 
sizes of order 10 11 cm, to very large structures, up to a 
few times 10 12 cm, i.e., 10-20 stellar radii. As discussed 
below, the largest events with a reliable size estimate 
(with small uncertainties) typically have L ~ 5 i?» . 

A number of intense flares on active coronal sources 
have been analyzed to date with the approach used 
here. These include a number of YSOs, both accret- 
ing and n on-accreting (YLW 15, HD 283572, LkHa,92, 
V773 Tau.lFavata et alJ200lD . the nearby ZAMS star AB 
Dor ijMaggio et alJl2000j) . a number of dMe flare stars 
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Fig. 2. — Top left panel: the characteristics spectrum of COUP 
1343, together with its best multi-temperature fit; top right: spec- 
trum collected during segment 2 (at the beginning of the flare rise) ; 
bottom left: spectrum collected during segment 4 (at the begin- 
ning of the flare decay); bottom right: spectrum collected during 
segment 9 (at the end of the flare decay). The flare spectra are 
shown together with the best-fit single-temperature model. 



iFavata et alJ l2000ri iFavata et all l2000al iReale et all 
120041) . and a number of active binary systems as well as 
Algol ijFavata fc Schmittlll999j) . In all cases, without ex- 
ception, the size of the flaring structure derived from the 
analysis is at most comparable to the stellar radius, and 
in most cases it's smaller. On the other hand the present 
analysis of intense flares on YSOs in the ONC results, 
in a number of cases, in very large loop sizes, extending 
to several stellar radii, showing them to be events of a 
different nature from the flares analyzed to date. 

The strongest evidence for large loop sizes in the 
COUP sample comes from well resolved long-lasting 
events with little sustained heating. In these cases, one is 
largely observing the undisturbed thermodynamic decay 
of the flaring structure, which, as discussed in Sect. |3J 
is directly related to the size of the flaring structure. 
One such example is COUP 1343, for which our analy- 
sis gives C = 1-95 ± 0.51, and F(() = 1.8, resulting in 
a loop size of ~ 1 x 10 12 cm. The presence of a small 
amount of residual heating is derived from the tail of the 
light curve, which decays more slowly, and indeed the 
detailed simulation of Sect. 14.11 shows that for the first 
50 or so ks the decay of both the emission measure and 
the temperature is well reproduced by a freely decaying 
loop of L = 1 X 10 12 cm. In fact, if we were to relax our 
conservative assumption on the peak temperature of the 
loop (with the best-fit temperature clipped to 100 MK) 
and use the nominal best-fit temperature, the resulting 
loop size would be even larger. 

Other similar cases of well determined flaring struc- 
ture sizes include COUP 669, 891, 971 and 1246. In all 
these cases -F(C) ^ 2.5, with a well constrained value of £, 
leading to a reliable determination of L from Eq.0 Both 
COUP 971 and 1246 have L~4x 10 11 cm, while COUP 
669 and 1343 have L~lx 10 12 cm. With L = 1.7 X 10 12 
cm, COUP 891 is the largest of the well constrained loop 
structures. In all these cases, the uncertainties for L re- 
sulting from the uncertainty in the determination of F (C) 
are of order 20-30%. Apart from COUP 971 (which, at 
L ~ 1.6 i?*, is likely to be a large "normal" coronal flare, 
with both loop footpoints anchored on the stellar pho- 
tosphere), the well determined loops have L ~ 4-5 i?*, 
which, as discussed in Sect. H3 is the typical corotation 
radius for a low-mass YSO (also likely to be the disk 
truncation radius), supporting the hypothesis that these 
loops may indeed be structures linking the star and the 
accretion disk. 

For the large events with a well constrained loop size, 
the peak densities are all in the range 1-8 x 10 10 cm~ 3 , 
and the equipartition magnetic fields (determined from 
Eq.EH a t the top of the loop) range from 40 to 250 G. In 
the assumption of a simple dipole geometry (very likely 
to dominate at large distance from the stellar surface) 
the photospheric field would be 



B 



(L - i?*) 3 x B c 



(11) 



Assuming an average equipartition field of 100 G 
and a typical loop length of 5 i?* , the correspond- 
ing field strength at the stellar surface would there- 
fore be of order 6 kG, at the top of the range of the 
field strengths determined in YSOs by Zeeman splitting 
l|.Tohns-Krull fc Valentill2005l) . 

In a few flares of Table Q] we obtain £ ^> 1.5 (e.g. 
COUP 332, 454, 848, 976), outside the range of validity of 
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Eq. [5j Probably, in these events the analyzed part of the 
decay is too short to cover a significant part of the decay 
path, and the slope is not yet well defined. Cases which 
result in long flaring structures in the presence of strong 
sustained heating are to be regarded with caution, as the 
strong heating dominates the decay, hiding the intrinsic 
decay of the cooling plasma and potentially introducing 
a larger uncertainty in the results. Also, the hyperbolic 
form of Eq. [S] will tend to include lengths much smaller 
than the 'nominal' results of Eq. within the final error 
bar for L. One such example is COUP 1410, for which 
C = 0.45 ± 0.12 and F(C) = 6.1. While the nominal 
loop length is L = 1.1 x 10 12 cm, the uncertainty range 
resulting from the uncertainty in £ is L = 10 n -2 x 10 12 
cm, with a factor of ten error bar at the lower end. 

A number of individual flaring events of particular in- 
terest are discussed in detail in Appendix [7| where the 
individual light curves and the time evolution of the spec- 
tral parameters are also shown. 

4.1. Detailed simulation of the flare on COUP 1343 

Although the method used here has been widely ap- 
plied in the literature, and its reliability tested for exam- 
ple in the case of the eclipsed flare observed by SAX o n 
Algol l|Schmitt fc Favatalfl999t iFavata fc Schmittlll999fl . 
as discussed in the above Section no previous analysis 
based on this approach has resulted in such large loop 
sizes (of order ~ 0.1 AU). It is thus legitimate to ask: 
1) whether such gigantic loops can be still described 
with standard loop models, closed magnetic structures 
where plasma is confined and moves and transports en- 
ergy along the magnetic field, and 2) whether the di- 
agnostics we are using (in particular the slope £ in the 
log T-log V EM plane) are still valid. To check on the ap- 
plicability of the method to such extreme regimes, and 
to validate its use on the present sample, we have mod- 
eled one of the well-resolved large flares from the COUP 
sample in detail. We have used a time-dependent hy- 
drodynamic model of plasma confined inside a coronal 
loop, as already done for other more regular stellar flares 
l|Reale et al.lll988tlReale'et~aIll200l . 

We have chosen to model the large event on COUP 
1343, in which the peak count rate is approximately 
30 x larger than the characteristic rate and which ap- 
pears, from the analysis performed using the approach 
described in Sect. 13.31 to have evidence for little sus- 
tained heating. The event lasts some 130 ks, and both 
the rise phase and most of the decay are not interrupted 
by gaps in the observations. Additionally, the decay in 
the log T-log V EM plane is well defined. The event and 
its analysis with the analysis method used here (as de- 
scribed in Sect. l3.3fl is shown in Fig.Q] 

To model the event on COUP 1343 in detail, we have 
assumed the plasma to be confined in a loop with con- 
stant cross-section and half-length L = 10 12 cm, sym- 
metric around the loop apex. The flare simulation is 
triggered by injecting a heating pulse in the loop which 
is initially at a temperature of ~ 20 MK (which is the 
dominant temperature in the characteristic spectrum as 
well as the temperature observed at the end of flare de- 
cay). We have chosen to deposit two heat pulses with 
a Gaussian spatial distribution of intensity 10 erg cm~ 3 
s _1 and width 10 10 cm (1/100 of the loop half-length) at 



50 100 150 
Time (ks) 



20 40 60 80 100120140 
Time (ks) 



3.0 - 
2.5 
g2.0 

| 15 
.2 1.0 
0.5 
0.0 _ 




20 40 60 80100120140 
Time (ks) 



-3.5 -3.0 -2.5 -2.0 
log(EM[10 57 cm- 3 ]) 



Fig. 3. — The results of the detailed hydrodynamic modeling of 
the large flare on COUP 1343. The top left panel shows the flare's 
light curve (as seen by ACIS), the top right panel the evolution of 
the emission measure, and the bottom left panel the evolution of 
the temperature. The bottom right panel shows the evolution of 
the flare in the log T— log \/EM plane. Continuous and dashed lines 
show the evolution of parameters from the hydrodynamic model, 
while the points with the vertical error bars are the values deter- 
mined from the time-resolved spectral analysis of COUP 1343. 

a distance of 2 x 10 10 cm from the footpoints, i.e. ver y 
close to the footpoints themselves l|Reale et al.l f20041 . 
After 20 ks the heat pulses are switched off completely, 
with no residual heating in the decay, as indicated by 
the density-temperature diagnostics for this specific flare. 
We have computed the evolution of the loop plasma by 
solving the time-dependent hydrodynamic equations of 
mass, momentum and energy conser vation for a com- 
pressible plasma confined in the loop l)Peres et al.lll982t 
iBetta et al.lll997j) . including the relevant physical effects 
such as the plasma thermal conduction and radiative 
losses (computed assuming solar metal abundances - at 
these high temperatures however the line losses are a mi- 
nor contributor, so that the results are robust against 
differences in the plasma metal abundance) . The gravity 
component along the loop is computed assuming a radius 
i?,* = 3 Rq and a surface gravity g* = 0.1 g®, typical 
of low-mass pre-main sequence objects. The results of 
the simulations are shown in Fig. |3J together with the 
observational results for COUP 1343. 

The computed evolution of the flare's parameters 
largely resembl es the evolut i on co mputed for other stel- 
lar flares (e.g. Rcale ct al. 2004), although on larger 
scales. The heat pulses make the loop plasma heat 
rapidly (~ 1 hour) to temperatures above 200 MK and 
the initially denser chromosphcric plasma near the loop's 
footpoints expands dynamically upwards at speeds above 
2000 km/s, to reach the loop apex on similar timescales 
(also about 1 hour). After this first impulsive phase, the 
temperature does not change much and the evaporation 
continues substantially but less dynamically. After the 
end of the heat pulse, the plasma begins to cool while 
still filling the loop. The density begins to decrease only 
two hours later. 

From the density and temperature distribution of the 
plasma along the model loop, we can synthesize the ex- 
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pectcd plasma X-ray spectrum filtered through the ACIS 
spectral response, and the relative spectral parameters as 
they would be determined through a spectral fit to the 
ACIS data. We have computed the emission in the coro- 
nal part of the loop, i.e. above the loop transition region, 
assuming a hydrogen column density AT(H) = 10 22 cm~ 2 
(i.e. the value obtained from the fit to the source's char- 
acteristic level spectrum). The dashed line in Fig. |3] (top 
left panel) shows the ACIS light curve resulting from 
the simulation, integrated above 1 keV, assuming a loop 
cross-section radius of 2 x 10 10 cm, i.e. (3 = 0.02. The 
smaller j3 than customary (/5 = 0.1 being assumed else- 
where) is necessary to provide a good fit to the flare's 
decay; the small f3 also results in footpoints still small 
relative to the photosphere. The spectra obtained from 
the loop modeling have been fit with 1-T model spectra, 
and the resulting evolution of the emission measure and 
of the temperature, as well as the log T-log V EM dia- 
gram are also shown in the figure. A visual comparison 
of the model results with the actual light curve and with 
the evolution of the measured spectral parameters indi- 
cates that the loop simulation results are in good quali- 
tative agreement with the data and that, therefore, the 
data are compatible with the flare occurring in a single 
giant flaring structure. The modeling could be further 
improved by changing some of the simulation's parame- 
ters; for example the light curve's agreement for t > 50 
ks could be improved by including a low resi dual heating 
activ e during the late decay (as done e.g. bv lReale et al.l 
120041 in their simulation of a large flare on Proxima Cen). 
However, such changes would not alter the large-scale 
features of the flaring loop, in particular its large size, 
and thus for the purpose of the present work they are 
not necessary. 

One additional diagnostic provided by detailed simula- 
tions such as the one discussed here is the rise time of the 
event. This is not available with the more general analy- 
sis employed for the complete sample of flares discussed 
in the present paper, which only uses the decay phase. 
In the case of COUP 1343 the rise phase is very steep, 
with a very peaked light curve. Not all impulsive heat- 
ing functions will produce such steep rise; we find that 
to produce the steep rise it is necessary to have strongly 
localized heating deposited very close to the loop foot- 
points. More diffuse heating, or heating localized near 
the loop apex would produce a smoother light curve. 

5. COMPLEX FLARING EVENTS 

While the present paper deals with all the flaring 
events in the COUP sample which can be model ed within 
the framework described bv lReale et al.l l|1997fl . a num- 
ber of events are present in the sample which do not 
show such regular behavior, and which cannot therefore 
be analyzed with the same approach. Perhaps the most 
dramatic one is the flare on COUP 450, which, coinci- 
dentally, was also observed at mm wavelengths. COUP 
45 is associated with a well known radio source (GMR- 
A, iGarav et al"lll987l) and is id entified with star HC 573 
from iHillenbr and fc Carpenter! l|2000ft . 

The X-ray event has been briefly reported by 
iCetman et atl {2003J) . and the COUP light curve for the 
event is shown in Fig. 0] As discussed in detail by 
IBower et all £2003), thc 

source also flared in the mm 
regime, although the mm flare started significantly later 
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Fig. 4. — Top panel: the COUP light curve of source 450. Bottom 
panel: evolution of the flare. See caption to Fig.^for details. The 
orange and green bars indicate the extent of the flaring intervals. 

(about two days after the start of the X-ray flare) . Thc 
mm flare lasted some 20 days, with erratic variability 
continuing for tens of days thereafter. The COUP ob- 
servation only covers the first day or so of the mm event 
(which has a much spottier time coverage) , and no obvi- 
ous correlation between the details of the two light curves 
is present. The counterpart to COUP 450 has been 
studied bv IBower et alJ l|2003j) who (using IR Keck spec- 
troscopy) identify it as a heavily absorbed K5 T Tauri 
star with average photospheric magnetic fields B ~ 2 — 3 
kG determined from Zeeman-splitted infrared lines. 

The COUP 450 event is not exceptional either in 
terms of its peak luminosity (with a rise of ~ 20 x in 
X-ray count rate at peak) or in terms of its duration 
(other events with durations of a few days are present 
in COUP). It is however one of the most complex events 
in terms of its light curve shape, and it lies at the ex- 
treme of complex flares, which lack a clear, well defined 
decay phase. Such events cannot therefore be analyzed 
with the method used here. In this specific case, the 
analysis is also hampered by the very strong degree of 
pile-up present in the Chandra observation, so that while 
the raw count rate ensures a 'high-resolution' view of the 
light curve, most events are lost to pile-up, and the resid- 
ual pile-up free events suffer from very limited statistics. 
Nevertheless it is possible to still resolve the light curve 
in a number of intervals, allowing us to monitor the evo- 
lution of the spectral parameters. 

The time evolution of the flare can be seen in the bot- 
tom panels of Fig. The very significant amount of 
pile-up in the ACIS observation results in a very small 
pile-up free count rate. In the light curve shown in the 
top panel of Fig. 0] the purple line shows the ML blocks 
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of the pile-up free events: although the raw count rate 
approaches 1 cts/sec, the pile- up free ML blocks have du- 
rations of tens of ks, which clearly subsample the event's 
evolution. Even so, however, the irregular evolution of 
the event is readily visible. No regular decay is present; 
the emission measure does not decrease, and the tem- 
perature shows an irregular behavior, with a strong in- 
crease in the next to last bin. Whether this is evidence 
of strong reheating within the same flaring structure or 
whether this is a separate structure undergoing flaring 
(perhaps induced by the first event) cannot be inferred 
from our data, specially because of the coarse time sam- 
pling caused by pile-up. 

The lack of a regular evolution for this event does not 
allow us to make (using the approach adopted in the 
present paper) any inference about the size and charac- 
teristics of the flaring region, and thus about the flaring 
mechanism and the underlying physics. The clear lack 
of a decay does point in this case to a complex heating 
function, with a highly time- variable release of energy in 
the flaring plasma, and - perhaps - to the involvement of 
different structures in the flaring, rather than to the sin- 
gle, confined plasma structure inferred from the analysis 
of the simpler events discussed in Sects. |21 and ^ 

6. DISCUSSION 

In a significant number of cases the size of the flar- 
ing structures derived here for intense flares in COUP 
sources is large (several stellar radii) with absolute sizes 
up to 10 12 cm. In a few of them, the evidence for large 
loop sizes is strong, because we could constrain negligi- 
ble heating to be present during the decay. Significant 
heating would have made the diagnosed lengths shorter 
and at the same time more uncertain. Given the high 
plasma temperatures at the peak of the flare (T > 100 
MK) the magnetic field necessary to confine the plasma 
typically is of a few hundred G. This scenario is very 
different from the one observed in more evolved stars, 
ZAMS or MS, in which the same type of analysis of flare 
decay as performed here has invariably yielded relatively 
compact magnetic structures L < R*, compatible with 
scaled-up version of the geometry observed on the Sun. 
Also the limited number of flar es previously studied with 
the same technique on YSOs l)Favata et al.l l2?KTjT) have 
shown evidence for compact structures. This is not in 
contradiction with the ev idence from the COU P sample, 
as the YSOs studied by iFavata et al.l ((2001) all are in 
the nearest SFRs (thus allowing studies of much less in- 
tense flares) and were the result of shorter observations 
(and therefore biased against the long lasting events well 
represented in COUP). 

Such long loops are unlikely to be anchored on the star 
alone, as the centrifugal force would rip them open and 
eject the plasma. The corotation radius is defined as the 
distance at which a test gas particle in keplerian orbit 
around the central star has the same angular velocity as 
the star, 

Rc = (GM./fi^) 1 / 3 (12) 

The effective gravitational potential (gravity plus cen- 
trifugal) takes the form of a vanishing gradient at the 
corotation radius, which is a saddle point, allowing only 
for magnetically dominated plasma to be present above 
and below the equatorial plane. The plasma favors mov- 



ing inward while inside of the corotation radius, but flow- 
ing outward while outside of the corotation radius. Un- 
less the magnetic field is strong enough to provide suf- 
ficient tension therefore enforcing corotation, the field 
itself would easily be opened up by plasma flowing down- 
hill in the e ffective gravitational po tential. 

Following iMestel fc SpruTtl <|1987|> - hereafter MS87 - 
and assuming the simplest stellar dipole field, a scaled-up 
version of the Sun for the stellar activity (with the mag- 
netic field strength proportional to the rotational velocity 
of the star, i.e. the p — 1 case of MS87, so that Lx oc B 2 ), 
Table 1 in MS87 gives a list of sizes of the largest closed 
loops that can possibly exist for a rotating solar-mass 
star under similar condition. By simply re-scaling with 
our stellar parameters, for loops of temperatures on the 
order of 100 MK or higher, the largest equatorial extent 
for closed magnetic loops (in the notation of MS87, C- 
cusp point, see their Fig. 1), will likely be in the range 
of 3-4 i?*, for a slow rotator (6-8 days, or a/a© ~ 4-5 
in the notation of MS87, where a is the rotation rate 
of the star), and high plasma temperature (£<j cx a and 
(Cd)o = 4, where Q is the ratio of magnetic to thermal 
energy in the corona) . Similar consi derations have re- 
cently been applied by Jardine (2004) in the context of 
'coronal stripping' in fast rotating, supersaturated stars. 

Our flaring sources have typical rotational periods 
(when they are known) of a few days, up to a week. 
For typical CTTSs periods of 6-8 days, the correspond- 
ing corotation radius is about 9 x 10 11 to 1.2 x 10 12 
cm, equivalent to ~ 5i?*, if i?* = 3Rq. Loop struc- 
tures larger than this would have difficulty anchoring on 
the star alone without very strong magnetic field ten- 
sion to counter balance the centrifugal force. What is 
then the magnetic geometry which can support coher- 
ent fields of hundreds of G over 10 12 cm? In the case 
of YSOs one obvious possibility is to have magnetic flux 
tubes connecting the stellar photosphere with the disk, 
at the co-rot ation radius or so mewhat outside it. As e.g. 
discussed bv lShu et alJ (|1997T> . the twisting of the mag- 
netic field lines induced by the differential (Keplerian) 
rotation of the inner disk rim and the photosphere will 
twist the flux tubes, resulting in longer loops and presum- 
ably in the stressed field configurations which drive the 
flaring. A good candidate among COUP flaring sources 
for star-disk flaring structures is COUP 1083, a moder- 
ate mass star with P rot = 5.9 d and with a flare implying 
L = 2.4 x 10 12 cm = 34i? Q , in the presence of mod- 
erate sustained heating, with an uncertainty range 1.9- 
3.0 x 10 12 cm. Even assuming a very young YSO, with 
a large radius R = 3i?0,the loop length is L ~ 10 i?* . 
Given the above consideration on stability of large mag- 
netic loops in rotating stars, it is unlikely that this type 
of loops can be anchored on the star alone. 

Magnetic flux tubes extending from the inner rim of 
the disk to the stellar photosphere provide a natural lo- 
cation for the intense flares found here on a number of 
sources. An additional open question is what heating 
mechanism can heat the plasma in these long loop to such 
high temperatures. While the general coronal heating 
mechanism is still not fully established, the most mod- 
ern m agn etohydrodynamic sim ulations (e.g. iPeter et alJ 
120041 and iSchriiver et alJl2004f) point to the heating be- 
ing due to shuffling of the loop's footpoint, a natural 
result of convective motions in the photosphere. The 
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photospheric convection is certainly present in accreting 
YSOs, and shearing of the disk would provide an equiv- 
alent, natural source of loop footpoint shuffling on the 
disk side. 

It should be stressed that the COUP flaring sources 
analyzed in the present paper present a large variety of 
flaring structures. While the long flaring loops which 
likely connect the star with the accretion disk are well 
represented in the sample, so are also the more compact 
structures which have very similar characteristics with 
the ones observed on more evolved stars, ZAMS or MS. 
Thus, it appears from the present analysis that two types 
of (flaring) coronal structures can coexist on YSOs: one 
which appears similar to the one present in older, more 
evolved stars (and thus likely to be a scaled-up version 
of the solar corona) and one made of very extended mag- 
netic structures, which is instead peculiar to YSOs, being 
dependent on the presence of disks on which to anchor 
one footpoint of the very long loops. While unfortu- 
nately none of the COUP sources has more than one 
flare sufficiently bright to be analyzed with the present 
technique, some sources (notably COUP 597 and COUP 
891, both discussed in Appendix show, in addition to 
the long events linked with long loops, some short events, 
too short to be analyzed, but likely originating in com- 
pact loops (given the small ti c ), so that it would appear 
that both types of coronae could exist, at the same time, 
in a given YS O. 

Wol k et all (|2005fl have presented an analysis of flare 
characteristics and frequency in "young suns" in the 
COUP database. Only two of our ob jects are in common 
with the sample of lWolk et aT]l)2005|) . namely COUP 223 
and 262. While no statistical conclusions can be drawn 
from only two objects, the flare on COUP 262 results in 
a long coronal structure (L > 3.6 see discussion in 
Sect. IA.3| . showing that these extended magnetic struc- 
tures were indee d present in the you ng Sun. 

Very recently, lLoinard et alJ l)2005|) have obtained ev- 
idence for the presence of long magnetic structures, in 
T Tau S: their radio VLBI observations show that, in 
addition to a compact source which they identify as the 
star's magnetosphere, a fainter streak of radio emission 
is present, extending to a distance of about 30 Rq from 
the star, which they state "may result from reconnection 
flares at the star-disk interface" . The size of the radio 
structure (> 10i?„, T Tau S having a radius R ~ 3i?©) 
is very similar to the loop lengths we have derived here 
from the analysis of the large X-ray flares, providing in- 
dependent support for the presence of such long, coherent 
magnetic structures in YSOs. 

6.1. The role of disks and of accretion 

To investigate the role of disks in the presence of long 
magnetic structures one would ideally like to have a 
clear separation of the sample in stars with disk and 
stars without. Unfortunately, the available diagnostics 
do not allow us to do this in an unambiguous way. The 
most complete investigat ion on the presence of d isks in 
the ONC is the work of IHillenbrand et all (|1998(K who 
have compared the results of a number of indicators for 
the presence of disks in the ONC, including IR excesses 
in a number of color bands and Can emission. Emis- 
sion in Ca II should indicate the presence of active accre- 
tion, and thus will miss "quiescent" disks. NIR excesses 



are often used as an indicato r of the presence of disks, 
and IHillenbrand et al.l <|1998(1 have quantitatively mod- 
eled the expected excesses in a num ber of bands , using 
models very similar to the ones of iMever et alJ l)1997fl 
and conv erging on similar conclu sions. As discussed in 
detail in [Hillenbrand et all l|1998f) the IR excess is ex- 
pected to be produced by the disk material being heated 
both by the stellar radiation and by the friction induced 
by accretion. In practice, its magnitude depends on a 
number of factors, such as the relative importance of ac- 
cretion vs. the size of the disk's inner hole, the relative 
contrast between the disk emission and the photosphere 
and the system inclination. Disks with low accretion rate 
and with a large inner hole, generating only a far-IR ex- 
cess could be difficult to see i n the NIR (see also the 
discussion m IMever et al.lll997D . 

Hill enbrand et al.l l)1998|i conclude that the excess 
A(I — K) (determined by taking into account the star's 
spectral type and thus effective temperature and intrin- 
sic colors and its reddening) is the most effective sin- 
gle indicator of the presence of disks; measurements in 
in the L band s would have be even b etter but were 
not available to H illenbrand et all l|1998j) except for few 
stars. However, while A(I — K) clearly is an effec- 
tive indicator on a statistical basis, its value for de- 
termining whether individual stars have a disk is much 
less clear. The d i stribu tions shown e.g. in Fig. 10 of 
Hille nbrand et alJ l)1998|) have strong tails to negative 
values (as far as A(I — K) ~ — 1) and a significant num- 
ber of stars with Can emission (thus, accreting systems) 
do not have A(J — K) exces ses. Even more significan tly, 
as discussed in Sect. 7.3 of IHillenbrand et al.l l)1998j) . a 
number of sources wit h disks visible in the H ST images 
("proplyds", from the lO'Dell fc Wondll996l sample) do 
not have IR excesses, clearly showing that some sources 
with disks escape detection through IR excesses. 

In the present paper, we have tried to analyze whether 
a relationship exists between the detection of long flaring 
structures and the presence of disks from the available 
indicators, bearing however in mind the above caveat. 
In Table ^ we report, for all the COUP sources ana- 
lyzed in our sample, the availa ble data from the compi- 
lation of IGetman et al.l 20051) - the original sou rce for 
most of the IR data is Hillenb rand et all l|1998j) - re- 
garding Can emission and excess in / — K colors. We 
also report whether the sources have an excess in the 
NIR JH KL bands from inspecti on of Fig. although 
(quoting Hill enbrand et alJ 119981 literally) "calculating 
disk fractions simply by counting the relative number 
of stars outside of and inside of reddening vectors in ob- 
served color-color diagrams is clearly naive" . 

Fig. plots the IR color-color diagrams for all COUP 
sourc es for which JHKL NIR photometry is available 
f from IGetman et al"]l2005j) . with the subsample studied 
in the present paper singled out. The slanted blue line 
in Fig. [5] is the ob servational mean lo cus of CTTSs in 
Taurus derived by IMever et al.l l)1997l) . All sources ly- 
ing to the right of the reddened photospheric locus in 
Fig. O (bounded by the green dashed lines) have NIR 
colors compatible (bearing in mind the above caveat) 
with the presence of disks similar to the ones present 
in Taurus YSOs. For objects whose NIR colors place 
them in the reddened photospheric region the situation 
is however much less clear as they could still have disks, 
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perhaps with low accretion and large interior holes, gen- 
erating negligible I R excesses. Indeed, as discussed by 
iMever et al.l l)1997j) - see their Fig. 5 - the most impor- 
tant parameter driving the NIR excess is the accretion 
rate, and for a given color excess there is a minimum 
accretion rate, thus disk heating rate, needed to gener- 
ate a NIR excess, independent of central hole size and of 
the disk inclination. Also, large holes (-fiWe 3> R*) will 
produce no NIR excess in the JHKL bands, with any 
excess shifted redward of 3 fim. 

From Fig.[S]it is evident that none of the COUP sources 
with strong flares has a significant H—K excess (with the 
exception of COUP 1608 and, marginally, COUP 597), 
while 6 sources (out of 13 for which L photometry is avail- 
able) show a K — L band excess. Of these 6 (COUPs 332, 
454, 597, 848, 976, 1343), 5 have no 77 - K band excess, 
while COUP 597 is the only one possibly displaying ex- 
cess es in both H — K and K — L. In the framework of 
the IMever et all <JI 
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models, COUP 597 is the only 
one whose JHKL colors would be compatible with the 
presence of an accreting disk with a moderate-size cen- 
tral hole, while the JHKL colors of COUPs 332, 454, 
848, 976 and 1343 (with K — L excess but no H - K 
excess) would be compatible with the presence of a cen- 
tral disk with a larger size central hole (few stellar radii) 
and low accretion rate. The remaining sources in our 
flaring sample have no measurable excess in these bands, 
which would be compatible cither with their being disk- 
less sources or their having disks characterized by a very 
low accretion rate and with large inner holes. 

An additional complexity comes from the fact that 
star-disk magnetic coupling requires a gaseous disk at 
the corotation radius, while NIR excesses are only sensi- 
tive to the dust content of the disk. But, given sufficient 
heating, the dust might be sublimated relatively near the 
star leaving behind the needed gas disk withou t an NIR 
excess. As discussed bv lMuzerolle et aD ()2003j) . for suf- 
ficiently large accretion rates, the inner rim lies beyond 
the corotation radius, so that pure gaseous disks must 
extend inside the dust rim. 

The situation is however clearly more complex, as ev- 
ident from the fact that some sources which fall, in the 
color-color diagrams of Fig. in the region of normal 
reddenened photospheres with no disks have significant 
A(7 — K) excesses. Unfortunately the overlap between 
the two samples is not very large (13 out of 32 sources 
have JHKL colors, 19 out of 32 have A(7 - K) de- 
terminations, and only 6 have both), but all 6 sources 
for which both indicators are available have a A (7 — K) 
excess, while 3 of them fall in the region of reddened 
normal photospheres in Fig. In the whole sample of 
1616 COUP sources for which optical data are available, 
472 have a A(J — K) determination; of these, 388 have 
A(7— K) > 0.02 (82%), compatible with an excess due to 
circumstellar material, while 71 have A(7 — K) < 0.02. 
In the present sample, A(7 — K) > 0.02 for 30 out of 
32 sources (94%) , a comparable number as the complete 
sample given the low statistics. 

Further evidence that disks may be present also when 
no excess in JHKL NIR colors is observed comes from 
the COUP data directly. iTsuiimoto et alJ 0005) have 
recently searched the COUP databases for sources show- 
ing significant Fe 6.4 keV fluorescent emission, of which 
they found 7 cases. In all cases, fluorescent emission 
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Fig. 5. — Top panel: the J - H vs. H - K NIR color-color 
diagram for all COUP sources for which the information is available 
(grey points). Individual sources analyzed in the present paper are 
marked as red points, with a label showing the COUP number. 
The purple lines are two isochrones, for 10 5 and 10 8 yr (nearly 
indistinguishable from each other), while the blue line represents 
the locus of CTTS. The two green dashed lines delimit the region of 
the plane which can be explained as reddened diskless stars; stars 
to the right of this region are showing a significant NIR excess 
indicative of the presence of a disk. Bottom panel: same plot for 
J — H vs. K — L colors. This diagram is a more sensitive probe of 
the presence of cooler disks, although many of them will peak in 
the far IR and thus will also not be visible in the L band. 

was detected during intense flares, a fact easily explained 
by the need for a sufficient number of hard (E > 7.11 
keV) photons to excit e the fluorescence. Two of the 
ITsuiimoto et alJ l)2005|) sources are also in our large flare 
sample, namely COUP 649 and 1040. The remaining 
5 fluorescing sources have too low statistics (either be- 
cause of the short duration of the flare or because of 
the low peak count rate) for our analysis. Flares are 
however not a necessary condition for fluorescent emis- 
sion, as shown by the detection of 6.4 ke V fluorescence 
in the p Oph YSO Elias 29 in quiescence bv lFavata et alJ 
(2005). At the same time, the large equivalent widths 
of the 6.4 keV line observed in the COUP sample (as 
well as in Elias 29) are compatible only (as discussed by 
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iFavata et all 120051 and iTsuiimoto et alJ 12005]) with the 
fluorescence originating in a centrally illuminated disk 
observed face on. A disk illuminated from above, or other 
topologies (such as fluorescence from the photosphere) 
would result in a much lower equivalent width, likely too 
low to be detected on COUP sources. Yet most of the 
ITsuiimoto et alJ l|2005j) sources (and in particular COUP 
649 and 1040), for which a disk is certainly present, given 
the observed X-ray fluorescence, show no JHKL NIR ex- 
cess, although some of them (e.g. COUP 649) do show 
a significant excess in A(J — K). Again, JHKL NIR 
excesses appear to be, in the ONC, a sufficient condition 
for a disk to be present, but not a necessary one, with 
A(I — K) providing a more sensitive indicator. 

Given the biases affecting the samples and the prob- 
lems affecting disk diagnostics, we refrain from formal 
correlation analysis on e.g. flaring loop size and the pres- 
ence of disk, limiting the analysis to qualitative consid- 
erations. Using A(J — K) as an indicator, the only two 
sources with A(I ~ K) < in the present sample (which 
we take to imply that disks are most likely not present) 
both have flares confined to compact magnetic struc- 
tures (COUP 7 with L = 0.1 R* and COUP 960 with 
L = 0.3 i?*), while both compact and large magnetic 
structures are present in the sources with A(J — K) > 0. 

The 6 sources with significant K — L excess comprise 
flares spanning both large loops and more compact struc- 
tures (comparable to the stellar size): large structures 
appear present on COUP 332 (although the determi- 
nation of loop size suffers from a large uncertainty), 
COUP 454, with L ~ 3 x 10 12 cm, COUP 848, with 
L ~ 2 x 10 12 cm, and COUP 1343, with L ~ 10 12 cm. 
On the other hand COUP 597 has L < 6 x 10 11 cm and 
COUP 976 has L < 8 x 10 11 cm. Large loops are also 
present in sources with no K — L excess, such as COUP 
669 (L ~ 10 12 ) and COUP 1083 (L ~ 2 x 10 12 cm), 
so that no clear correlation appears present between the 
presence of large loops and K — L excess. With one sin- 
gle exception (COUP 1608) none of the sources in the 
present sample show significant H — K excesses. 

Only one strongly accreting source (COUP 141, with 
Can equivalent width —17.8 A) is present in our sam- 
ple of large flares; two sources (COUP 1114 and COUP 
1608) show very modest emission (and presumably ac- 
cretion rates), and all other sources show (when avail- 
able) Can in absorption or filled in and thus little or 
no ongoing accretion. The statistical significance of this 
is difficult to assess given the bia ses of our sample (dis- 
cussed in Sect. |2Jl, and the fact ijPreibisch et al.ll2005(l 
that accreting sources in the COUP sample are statis- 
tically intrinsically fainter than the non-accreting sam- 
ple (as well as more absorbed), introducing a small bias 
against strong flares from accreting sources being present 
in our sample. Interestingly, COUP 141 is one of the 
cases of compact loops, with L < i?*, so that in this 
case no disk-photosphere magnetic structures need to be 
postulated, and the flare is likely to be a normal coronal 
event, similar to the ones observed in more evolved stars. 
COUP 1114 and COUP 1608 both only result in upper 
limits to the size of the flaring structure. 

In the complete COUP sample, out of 1616 sources 
for which optical data are available, Can data are avail- 
able for 537 sources. Of these, 198 have Call in absorp- 



tion, 189 have Can in emission and 150 have filled lines 
(equivalent wi dth of the Call lines reported as "0.0" in 
iGetman et aT]|2005|l . Among the sources with Can de- 
termination the fraction of sources in emission is then 
35%, while in the flaring sample discussed in the present 
paper only 3 sources out of 20 for which the data are 
available have Can in emission, i.e. 15%. Given the small 
numbers, and the biases present in the various samples, 
it's difficult to assess the significance of the difference. If 
the lack of long flaring structures in accreting YSOs is 
however real, given that long magnetic structures con- 
necting the disk to the star are certainly present, it im- 
plies that active accretion inhibits the formation of strong 
flares in the extended magnetic structures connecting the 
disk to the star. Long and intense flares (and thus very 
hot plasma) appear to be present in non-accreting struc- 
tures connecting the disk to the star. 

7. CONCLUSIONS 

In the complete COUP sample of YSOs in the ONC, 32 
flares have sufficient statistics (in terms of flare duration 
and photon count rate) and regular evolution to grant a 
detailed analysis of their decay. The type of flare decay 
analysis employed here allow us to derive the size of the 
flaring loops in the presence of heating extending into the 
decay phase, and, under simple assumptions, estimates 
of the plasma density and confining magnetic field. 

The most notable result from our analysis is the strong 
evidence for very large flaring structures in these stars. 
The magnetic structures confining the plasma in a num- 
ber of cases are much larger than the stars themselves. 
Among the events analyzed, a large fraction have very 
high peak temperatures (T > 100 MK), and some are 
very long-lasting, with the longest flares extending to up 
to a week in duration. The peak X-ray luminosity for the 
most intense flare reaches 8 x 10 32 erg s _1 , comparable 
to other very bright flares observed in YSOs. 

Our sample is limited to the brightest 1% of COUP 
flares, and the results may not be representative of the 
average magnetic reconnection event in Orion stars. Nev- 
ertheless, the present results show that very extended 
magnetic structures confining hot flaring plasma exist in 
YSOs. Structures of comparable sizes have never been 
seen in more evolved stars, and, given the short rotational 
periods of many of the flaring COUP stars, would not be 
stable if anchored onto the photosphere with both foot- 
points. These long magnetic structures are in the present 
paper interpreted as linking the stellar photosphere with 
the inner rim of the circumstellar disk. 

The available indicators of the presence of disks do 
not allow for a formal analysis of whether large flaring 
loops are linked to the presence of disks; in particular 
stars without NIR excesses are not necessarily diskless, 
as they may e.g. have disks with little dust. However in a 
limited number of sources the NIR data clearly point to 
disks being present, either as K — L excess in color-color 
diagrams or as strong excesses above the de-reddened 
photospheric colors (A(I — K) > 0.4). Both small flaring 
structures (likely anchored on the photosphere only) and 
long flaring structures, with sizes which would extend 
well into the disk (assuming the disk to be truncated at 
the co- rotation radius) are present in these stars. Both 
types of structures (loops anchored on the star only and 
loops connecting the star to the disk) are then perhaps 
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likely to coexist in YSOs with disks. 

For the two sources in our sample without any A(I—K) 
excess, and for the star with a strong Can accretion 
signature, the flaring structures are relatively compact. 
This suggests that inner disks are needed for reconnec- 
tion in star-disk magnetic loops and that heating to X- 
ray temperatures is inhibited by mass loading associated 
with YSO accretion. More statistics are clearly needed 
to corroborate the present conclusions, and hopefully fu- 
ture, already approved long X-ray observations of other 
star-forming regions will be able to supply the needed 
additional evidence. 

The data presented in the present paper thus consti- 
tute the first direct evidence for the presence of very long 
magnetic structures in YSOs, and, by inference, of mag- 
netic structures linking the stellar photosphere with the 
circumstellar disk. Such structures are postulated by 
magnetospheric models of YSO accretion, but had not 
been directly detected to date. 
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APPENDIX A: NOTES ON SELECTED INDIVIDUAL 
EVENTS 

The present Appendix presents a detailed discussion 
(together with the relevant light curves and diagrams 
from the analysis) for a number of individually selected 
COUP sources. Its purpose is both to discuss in detail 
some peculiar or particularly representative event and 
to show the variety of flare shapes and characteristics 
present among the ONC YSOs. 

A.l. COUP 28 

COUP 28 is a 0.53 M , 2.3 R® star of spectral type 
MO, showing a NIR excess (A(I — K) = 0.30) and no ev- 
idence for active accretion. Its rotational period is 4.4 d. 
It displays one of the longest flares (Fig. EJ) of the COUP 
sample, with the total duration of the event extending 
over more than a week. Unfortunately an observation 
gap has blocked much of the flare rise, which however 
appears to be slow and extending over more than one 
day. The peak seems to be observed, and this allows a 
reliable analysis of the decay phase. The sources is at a 
large off-axis angle, so that no pile-up is present even if 
the flare at peak reaches a high count rate. The flare's 
peak count rate is almost 100 x the source characteristic 
level. The peak temperature, at ~ 80 MK, is a moderate 
one for COUP flaring sources (although it is very high by 
the standard of flares observed in older stars), and is still 
well determined on ACIS spectra. The log T vs. log V EM 
diagram shows, as typical in both solar and stellar flares, 
the temperature peaking well before the EM. The ex- 
cellent stastistics allows us to follow the flare's evolution 
in detail, and the temperature decay shows evidence for 
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Fig. 6. — Top panel: the light curve of COUP 28. Bottom panels: 
evolution of the flare (see Fig. for details). 

reheating (in block 9, Fig. |HJ), as also reflected in the 
change in slope of the light curve decay, which is not 
a simple exponential but rather shows evidence for two 
different time scales. In this case, fitting an average de- 
cay to both the light curve and the slope in the logT 
vs. log V EM diagram still provides a good estimate o f 
the parameters of the flaring region l)Reale et alJ l2004h 
At L = 5.5 x 10 11 cm = 1.9 i?* the flare is a relatively 
compact one among the ones observed in COUP. 

A.2. COUP 43 

COUP 43 is a known SB2 binary ; the mass and ra- 
dius computed by Hillenbrand (1997) for the unresolved 
source are 0.40 M and 2.9 i?©. Its spectral type is Ml, 
and it shows no evidence of active accretion; an I — K 
excess (A(7 - K) = 0.50) is present. The flare (Fig. 0) 
is a double event, with the second flare beginning while 
the first one is still in its decay. The rise phase of the 
first event is lost in an observation gap, and therefore 
only the second event, with a duration of about two 
days, and which shows an impulsive rise, has been an- 
alyzed. While the limited statistics only allow us to sub- 
divide the decay in three intervals, the steep £ implies 
very limited sustained heating, and thus a long structure 
(L = 1.1 x 10 12 cm = 5.5 R*) driven by the slow decay 
(ti c = 62 ks). The two flares are sufficiently well sepa- 
rated, so that the first (incomplete) event does not sig- 
nificantly affect the spectral analysis of the second event. 
However, the lack of coverage of the latest phases of the 
decay imply caution in the interpretation of the analysis 
results, as a flattening of the temperature decay could 
imply a final shallower £ and thus a smaller loop size. 

A.3. COUP 262 
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Fig. 7.— Flare evolution of COUP 43. 

COUP 262, with 1.1 M and 1.6 R e is a "young Sun" , 
and thus of particular interest for flare studies. In the 
co ntext of COUP, it s flaring activity is also discussed 
bv lWolk etaLl (|2005ft . With spectral type K5, it shows 
no evidence for active accretion and has a rather strong 
I — K excess (A(7 — K) = 2.2), pointing to a relatively 
massive but inactive disk being present (although no 
H — K excess is present). Although the flare (Fig. [HJl 
has a very sharp peak, it decays slowly {t\ c = 120 ks), 
lasting for a few days. The event represents a relatively 
modest increase over the characteristic level of a factor 
of 4, but the shallow C, implies a high level of sustained 
heating during the decay, so that (also given the high 
temperature, T > 100 MK) the resulting flaring loop is 
very large, L = 2.0 x 10 12 cm = 18 although the sig- 
nificant uncertainty in £ implies a lower confidence range 
of L = 3.6 Nevertheless, even at the lower end of 
the confidence range, the flaring region is larger than 
normally found in active stars, providing evidence that 
large loops, likely connecting the star to the disk, must 
have existed for the young Sun. 

A.4. COUP 597 

COUP 597 is a somewhat more massive star than 
COUP 262, with M = 1.5 M Q and R = 2.0 R Q . The 
strong Ca II absorption shows it not to be actively accret- 
ing. The flare (Fig.|5J) appears to be the superposition of 
two events, a less intense and more long-lasting one and 
a more intense one of which only the very fast decay is 
visible, as the rest falls into an observation gap. As the 
analysis relies on the delay of the second, slower flare, 
the first, fast event should have no effect on the results. 
The peak temperature for the longer event (determined 
during the rise phase) may be affected at some level from 



Fig. 8.— Flare evolution of COUP 262. 

the presence of the first flare. However, given the very 
moderate peak temperature, and the weak dependance 
(as \/T) of the loop size on peak temperature, this will 
not have a strong effect on the results. 

Of interest here is the slow rise of the longer flare, 
which appears to take place over a day or so. The loop 
size resulting from the analysis is L = 1.6 R*, thus a 
moderate size loop. The size is largely driven by the slow 
temperature decay and resulting shallow £ (with -F'(C) = 
9.7, implying strong sustained heating). The slow rise 
can therefore imply either a very slowly rising heating 
function, or perhaps (using the results of the modeling 
of COUP 1343, Sect. I4.1JI a heating source distributed 
along the loop (or even concentrated near the loop top) 
rather than localized at the loop footpoints. 

A.5. COUP 649 

COUP 649, is a typical low-mass member of the ONC, 
with 0.4 M Q and 2.2 i? Q . It shows no evidence for active 
accretion and has an / — K excess (A(I — K) = 0.4). 

The flaring event has a very peculiar evolution 
(Fig. llOjl . with an initial very slow rise, lasting for about 
a day, followed by a plateau at about 10 x the character- 
istic rate; after another half day, an impulsive peak in- 
creases the count rate by another factor of two, and then 
the decay begins. Notwithstanding the oddly shaped 
light curve, the event shows some of the characteristics 
of a magnetically confined flaring structure, in particular 
the temperature which peaks while the count rate is still 
increasing, and has already started to decrease on the 
plateau. Somewhat peculiarly, the secondary sharp peak 
preceding the decay does not result in an increase in the 
temperature. 

The limited statistics does not allow a detailed analy- 
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Fig. 9.— Flare evolution of COUP 597. 

sis, with only two points on the flare's decay. The peak 
temperature (T = 80 MK) is moderate by the standard 
of the other large flares analyzed here, and the nominal 
length of the flaring structure is L = 6.4 x 10 11 cm = 
4.2 i?*, although with a large uncertainty driven by the 
uncertainty on £. 

A.6. COUP 669 

COUP 669 is a relatively massive object, with 1.5 M© 
and 2.6 R Q , and an I - K excess (A(7 - K) = 0.4). The 
flare f Fig. II If) is a well defined impulsive event, for which 
the rise phase has however been cut by an observing gap. 
The peak temperature, at T — 79 MK is relatively mod- 
est, but the fast decay of the temperature implies limited 
sustained heating, resulting in a flaring loop with a rea- 
sonably well determined size, L = 9.2 x 10 11 cm = 5.1 i?*, 
typical of the large flaring loops in the COUP sample. 

A.7. COUP 752 

Another typical low-mass object, COUP 752, has 
0.5 Mq and 1.7 i?©, no evidence for ongoing accretion, 
and anl-K excess (A(J - K ) = 0.6). COUP 752 un- 
dergoes a very intense flare (Fig. I12|) . with a peak count 
rate more than 100 x the characteristic one. The event 
lasts over 5 days, and the rise phase is long (half a day) 
and well observed. The excellent statistics allow for a de- 
tailed analysis of the event, although the initial part of 
the decay unfortunately falls in an observing gap. The 
flare's peak temperature slightly exceeds 100 MK, still 
falling in the domain in which the ACIS detector gives 
reliable temperature determinations. 

Notwithstanding the very large increase in count rate, 
the event follows well the evolution in the log T-log V EM 
plane for a flare confined in a single loop. The temper a- 



FlG. 10.— Flare evolution of COUP 649. 
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Fig. 12.— Flare evolution of COUP 752. 

ture peaks (block 1) when the emission measure has just 
started to rise (and still is only about one tenth of the 
peak value), and then decreases to ~ 50 MK when the 
emission measure has peaked. The decay does not follow 
a simple exponential law; rather, there is clear evidence 
for sustained heating, with well defined reheating events. 
In particular, in block 11 (immediately after the first ob- 
serving gap) the temperature increases again briefly, and 
the decay of both the light curve and the emission mea- 
sure changes slope, clear evidence for reheating of the 
plasma. 

Our analysis shows a very shallow £, indicative of the 
presence of strong ongoing heating during the flare decay, 
as it is also evident from the increase in plasma tempera- 
ture. As a consequence, only an upper limit to the flaring 
structure's size can be derived, at L < 7.4 x 10 10 cm, with 
a best-fit value L = 6.5 x 10 10 cm = 5.6 i?*. The result- 
ing minimum confining magnetic field, at B > 250 G is 
relatively strong. 

A.8. COUP 891 

COUP 891 is the most massive object in the COUP 
flaring sample, with 2.4 M@, 4.9 Rq, a strong I — K ex- 
cess (A(J — K) = 1.10), no evidence for accretion and 
significant obscuration (Ay = 8.0 mag). The X-ray light 
curve suffers from significant pile-up, and therefore we 
have used, for our analysis, the list of photons extracted 
from an annular regi on free from pile-up, as described by 
iCetman et al.l l)2005fl . This strongly reduces the number 
of photons available for a spectral analysis, so that only 
moderate statistics are available for this event. 

Two flares are present in the COUP light curve 
(Fig. I13|) : however the first one is too short and can- 
not be analyzed with the approach adopted here. The 



Fig. 13.— Flare evolution of COUP 891. 

second event is a well defined impulsive event, with a 
fast rise and a regular, exponential decay, well traced for 
more than four days. At the peak, the flare count rate is 
~ 50 x larger than in quiescence. The peak temperature 
is slighly above 100 MK, but still within the range of the 
ACIS instrument. The fast decay of temperature results 
in a fairly steep £, and thus in little sustained heating. 
The long decay time (t\ c = 73 ks) therefore results in a 
long flaring structure, L = 1.7 x 10 12 cm = 5.1 i?*. 

A.9. COUP 971 

COUP 971 is a low mass (0.7 Mq) star with R = 
3.3 Rq. The Can in absorption points to a non-accreting 
system. The flaring event (Fig.ll4|l is a typical example of 
a short impulsive flare. The decay in the logT-log \JEM 
plane is steep and regular, pointing to little if any residual 
heating. The result is a moderate size loop (L = 1.6 R*), 
similar to the ones resulting from the analysis of intense 
flares in more evolved stars; this is likely an example of 
an event confined in a structure anchored to the stellar 
photosphere only. 

A. 10. COUP 1083 

Little optical information is available for COUP 1083, 
and in particular no mass and radius estimate is known. 
A rotational period has however been derived (P = 5.9 d) 
for this star. Its spectral type has been estimated at M0, 
showing that the object must have a moderate mass. 

The flare present in the COUP 1083 lightcurve has a 
peculiar shape, well visible in Fig. 1151 with a very slow 
rise, lasting more than one day, and a decay lasting about 
3 days. The peak of the flare unfortunately falls in an 
observing gap, and it thus not visible. We have analyzed 
the event assuming that the observed maximum is the 
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Fig. 14.— Flare evolution of COUP 971. 

actual peak, but in practice the peak might be higher, 
although this would not substantially affect our conclu- 
sions. 

Notwithstanding the peculiar light curve, the evolution 
of the event follows the one theoretically expected for a 
confined flare, with the temperature peaking prior to the 
EM and a regular decay in the log T-log V EM plane. 
A second event may be starting in the final phase of 
the decay of the flare, on day 12, as shown both by the 
increasing count rate and by the increasing temperature. 
The second flare is cut by the observation gap and thus 
nothing can be derived for it. We have not included the 
last ML block in the analysis. 

Although the statistics are limited (and thus the error 
bars large) , the derived steep £ is compatible with limited 
sustained heating, so that a long structure, L = 2.4 x 10 12 
cm, is derived from the analysis due to the slow flare 
decay (t\ c = 125 ks). The large size is also compatible 
with the very slow rise time, as discussed for the detailed 
model of COUP 1343 fSect.Ol). 

A.ll. COUP 1246 

COUP 1246 is one of the lowest-mass stars in the 
present sample, with 0.2 M Q and 1.6 Rq. It shows no ev- 
idence for active accretion, but it has a sizable K-band 
IR excess (A(7 — K) = 0.8) and is one of the sources 
showing an K — L excess but not a H — K band one; its 
rotational period has been determined at P = 5.2 d. The 
flare fFig. H6J) is a very well defined impulsive event, with 
a fast rise and a peak count rate ~ 100 x the characteris- 
tic rate. Also in this case the peak observed temperature 
is somewhat above 100 MK, but still within the range 
of reliability of ACIS. The decay in the log T-log \[EM 
plane is well defined, and closely follows the one expected 



Fig. 15.— Flare evolution of COUP 1083. 

from the models of confined flares, allowing to analyze 
the event in some detail. The light curve plotted in the 
top panel of Fig. ^] shows clear evidence for a second, 
smaller flare developing on the tail of the first one, re- 
sulting in an increase in temperature at the flare end. 
Statistics are insufficient to analyze the second flare at 
any level of detail, and therefore the last block has not 
been included in the analysis. The resulting size of the 
flaring region is L = 4.0 x 10 10 cm = 3.8 i?*, i.e. a typical 
event for the COUP sample. 

A.12. COUP 1343 

Very little is known about the optical counterpart to 
COUP 1343 (the analysis for the flare is shown in Fig.^l, 
with no mass or radius estimates known, although a rota- 
tional period measurement (P = 8.7 d) is available. The 
flare is a very well defined event, with a clean impulsive 
rise and an undisturbed decay. The peak flare luminosity 
is > 10 x larger than the characteristic source luminosity, 
allowing a clean analysis of the flare spectral character- 
istics. In the analysis we have ignored the last segment 
of the flare (the one extending into day 13), for which 
the temperature has too large an error (also because of 
the data gap) . The peak temperature of the flare is very 
high, with a nominal best-fit temperature of ~ 400 MK. 
As discussed above such high temperatures cannot be re- 
liably determined with ACIS, and therefore the resulting 
flare size (based on the assumption T b s = 100 MK) is 
to be considered as a lower limit. 

The rapid temperature decrease implies that little sus- 
tained heating is present, so that the observed decay 
is close to the thermodynamic decay of the loop. The 
smooth decay also argues for the lack of significant re- 
heating. Therefore, with t 1c = 39 ks and T obs > 100 MK, 
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Fig. 16.— Flare evolution of COUP 1246. 



§ 40 - 

s 



10 R : 0.51 ' Ets: B154 



0.01 0.1 
VEM [10" cm- 3 ] 




5x10 s 

s. since obs. start 



Fig. 17.— Flare evolution of COUP 1410. 



the analysis results in L ~ 1 x 10 12 cm, or L ~ 14 Rq. 
Even if COUP 1343 were a relatively massive star, with 
a radius of a few Rq , the flaring structure would still be 
several times the size of the star. The resulting magnetic 
field, while moderate in itself (150 G), must be coherently 
organized over some 0.1 AU. 

As described in Sect. 14.11 the evolution of the same 
event has also been simulated in detail, confirming that 
the temperature and emission measure evolution are the 
result of a flare confined in a large loop of the same size 
of the one found by the decay analysis discussed here, 
with impulsive heating. 

A.13. COUP 1410 

COUP 1410 appears more evolved than the rest of the 
sample, with an estimated age t = 3.6 x 10 7 yr, and a con- 
sequently small radius (0.51 Rq) for its mass (0.36 M©) 
It shows no evidence for accretion from Can emission 
but it has a significant if-band excess (A(7 — K) — 2.3) 
The rotational period has been measured at P = 6.76 d 

The flare (Fig. I17|) has a well defined impulsive rise 
followed by a very long decay extending for over a week 
with Tic = 150 ks. Even though significant heating is 
present (as shown by the shallow Q , the resulting size of 
the flaring structure is L = 1.1 x 10 12 cm = 55 very 
long with respect to the (small) stellar radius. However, 
the uncertainty on £ results in a significant uncertainty 
on L, which has a lower range of L = 1 x 10 11 cm = 5 
the same large size as found for the best determined large 
events in the COUP sample, making this event a good 
candidate for a star-disk connecting structure. 
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Fig. 18.— Flare evolution of COUP 1568. 
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COUP 1568 is another relatively massive object in the 
flaring sample, with 2.6 Mq and 4.0 Rq. It has a small 
IR excess (A(J - K) = 0.2). 

The flare (Fig. I18|l is relatively short, lasting just one 
day, but it is intense, peaking at over 100 x the character- 
istic rate, with excellent statistics due to the source's off- 
axis angle which prevents pile-up. The impulsive phase is 
very sharp, and the decay almost a perfect exponential. 
The peak temperature is very high: the best-fit values 
for the ACIS spectra at peak are T ~ 500 MK, however 
with a very large uncertainty which always include values 
as low as 100 MK, due to the ACIS spectral response. 

Once in the decay phase, the temperature decays very 



slowly, resulting in a very shallow £. Therefore only 
an upper limit to size of the flaring structure can be 
obtained, L < 1.5 x 10 11 cm, with a best fit value 
L = 3.7 x 10 10 cm = 0.4 R*. Unlikely many of the 
large COUP flares studied in this paper, this event is 
very likely to be compact, with size smaller than the 
star itself. A small size is also implied by the very fast 
rise phase (which constrains the time scale with which 
the chromospheric evaporation will fill the flaring loop). 
Given the compact size, together with the large peak 
emission measure and high temperature, the peak flare 
density will likely be high (n e — 1.2 x 10 12 cm -3 ), as 
must the confining magnetic field (B — 3.5 kG) 
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